The interior spacetimes of stars in Palatini f{R) gravity 
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We study the interior spacetimes of stars in the Palatini formalism of f{R) gravity and derive 
a generalized Tolman-Oppenheimer-Volkoff and mass equation for a static, spherically symmetric 
star. We show that matching the interior solution with the exterior Schwarzschild-De Sitter solution 
in general gives a relation between the gravitational mass and the density profile of a star, which 
is different from the one in General Relativity. These modifications become neglible in models for 
which 5F{R) = df /dR — 1 is a decreasing function of R however. As a result, both Solar System 
constraints and stellar dynamics are perfectly consistent with f{R) — R — /R. 

PACS numbers: 98.80.-k, 95.35.+d, 04.50.+h 



I. INTRODUCTION 

Current observations indicate that the expansion rate 
of the universe is accelerating [l|, 0] , which in Einstein's 
theory of General Relativity (GR) implies the existence 
of a non-zero vacuum energy density, the cosmological 
constant. Predicting a finite, small vacuum energy is a 
very serious problem from the quantum field theory point 
of view, in addition to which the cosmological constant 
model requires huge amount of fine-tuning to be able to 
explain the observed acceleration at present. One way 
to alleviate these problems is to assume a zero vacuum 
energy density (due to some yet unknown reason) and 
instead modify the gravitational action so that dynamics 
may give rise to the observed acceleration. One popular 
class of such models consists of generalizing the gravi- 
tational action to also contain non-linear interactions in 
the Ricci scalar R: 



S = 



1 



167rG 



-gCiAg^iy) , (1) 



where >C„i is the matter Lagrangian. Setting f{R) — R — 
2A corresponds to the canonical Einstein-Hilbert action. 
It has been shown that the addition of a 1/R term to 
the Einstein-Hilbert action can give rise to the observed 
acceleration at present [1, . 

In this approach one naturally also has to check if the 
modified gravity theory is compatible with other obser- 
vations. Indeed, as was shown by Chiba using the 
equivalence between /(i?) gravity, understood as a pure 
metric theory, and an w = Brans-Dicke model, a 1/R 
modification to the Einstein-Hilbert action is not com- 
patible with the post-Newtonian (PPN) constraints from 
Solar System measurements. This conclusion was con- 
firmed, after a recent debate in the literature [1, by 
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Erickcek, Smith & Kamionkowski [qI], who presented an 
explicit solution of the field equations arising from the 
action llj when considered as a pure metric theory. The 
conclusion of Erickcek et al. has since been shown to also 
apply for more general forms of f{R) 10}. 

In the metric formalism one assumes that the affine 
connection of the spacetime manifold is the Levi-Cevita 
connection, that is, the Christoffel symbol of the second 
kind {P,}: 



1 



(2) 



This choice corresponds to making sure that the met- 
ric geodesies, corresponding to the shortest (or strictly 
speaking extremal) paths between two points on the man- 
ifold, coincide with the the affine geodesies, i.e. the curves 
for which the tangent is always parallel transported to it- 
self. 

However, the restriction ([2|) is by no means obligatory, 
and more general spacetimes can be studied if the the 
affine connection is allowed to be a free variable inde- 
pendent of the metric. Not surprisingly, the solutions 
to the equations of motion arising from the action ([T]) 
in this larger space of manifolds, are in general different 
from those obtained in the restricted set of spacetimes 
corresponding to ([2]) . This extension of the configuration 
space is what is usually called the "Palatini-formalism" , 
and it is the point of view to f{R) gravity taken in this 
paper. We shall refer to this extension of the possible 
spacetimes as "metric-affine" configuration space, in con- 
trast with the "g-metric" configuration space correspond- 
ing to (21). For earlier work on gravity in the Palatini 
formalism, sometimes also referred to as the first order 
formahsm, see fsl. [TTI. [T2| . 

Varying the action ^ with respect to both the metric 
g^iy and the affine connection V^^, results in the equations 
of motion for Palatini f{R) gravity: 



F{R)R^,--f{R)g^, = 8^GT^,, (3) 
^p{^lF{R)gPn = 0, (4) 
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where F = df/dR, R = g^"'R^,^(T), T^^ is the energy- 
momentum tensor of and is the covariant deriva- 
tive with respect to the afhne connection F^^. Note 
that (H]) reduces to the g-metric compatibihty equation if 
F{R) = const. In this case, and only in this case, does the 
extremal solution for the full metric-afhne configuration 
space happen to reside within the subspace of (7-metric 
manifolds. The free affine connection then reduces to the 
usual Levi-Civita connection ([2]) and the f{R) model col- 
lapses to the usual Einstein-Hilbert action, possibly with 
a cosmological constant. 

As was emphasized by Erickcek et al. when studying 
metric f(R) gravity in Ref. ^9||, a unique exterior solu- 
tion for a stellar object is found by matching it with an 
interior solution in the presence of the matter sources. 
Indeed, the failure to find a metric interior solution that 
could be matched to an asymptotic Schwarzschild-De Sit- 
ter spacetime is the reason why most f{R) theories of 
gravity are ruled out in the metric formalism p^ . Mo- 
tivated by this observation we will here study the com- 
plete spacetimes of spherically symmetric static stars in 
the Palatini formalism. The hope is of course that new 
solutions can be found beyond the subspace of g-metric 
manifolds, where the interior solution can be matched 
to an exterior metric. This indeed turns out to be the 
case [l3| . 

We will in this paper match the spherically symmetric 
vacuum solution to the stellar interior in the Newtonian 
limit. In particular, section II will present an explicit so- 
lution for the case of a constant density. In section III, we 
will derive the exact Tolman-Oppenheimer-Volkoff and 
mass equation for f{R) gravity in the Palatini formal- 
ism, from which the equilibrium configuration of a rela- 
tivistic star can be computed once the equation of state 
is given. Unlike the case in GR, the mass equation can 
not be expressed as a single integral, but is instead cast 
in the form of a differential equation. We find that in 
the Palatini formalism, the stellar interior can always be 
matched to an exterior Schwarzschild-De Sitter solution, 
whereby these models easily pass the Solar System PPN- 
constraints However, we also find that Palatini f{R) 
models lead to a different relation between the appar- 
ent mass of the star and the interior density profile than 
does General Relativity. As a result, it will in principle 
be possible to constrain Palatini f{R) models based on 
astrophysical considerations. However, it turns out that 
for any model for which 5F{R) = F{R) — 1 is a decreas- 
ing function of _R, these corrections are neglible and the 
stellar interior becomes very similar to the usual solution 
in GR. In particular, the model with 5f{R) ^ 1/i? is 
perfectly consistent with both the existing Solar System 
constraints and the stellar astrophysics. 



II. CONSTANT DENSITY SOLUTIONS 

Let us first consider the solutions to the field equations 
when the density p is constant and pressure p can 



be neglected. Taking the trace of ([3]) one finds a purely 
algebraic equation relating R and T: 



F{R)R - 2f{R) = SttGT 



(5) 



Obviously, in the case of a constant T(p), the Ricci scalar 
is also a constant and the full field equations ([3J3]) are 
equivalent to the regular Einstein equations with a scaled, 
p dependent gravitational constant G/Fp-. 



(6) 



where G^^ = R^v—^Rpg^.i' and the effective cosmological 
constant Ap is given by 



A. = - f i?. - ^ 



(7) 



Since ([4]) now implies that the affine connection is given 
by ([2]), the Einstein tensor G^^ can be computed in the 
usual way from 5^^. 

We now want to study the general, spherically sym- 
metric solutions for the metric g^^, which can be written 
in the form: 



5^(*.-)dt2 + eS(*,-)dr2-i-r2dr!2. (8) 



Let us first consider the exterior solution. Because the 
relation between R and T is purely algebraic, the interior 
solution will not affect the vacuum value of i? = i?o- 



Furthermore, since in vacuum T,, 



0, the problem is 



reduced to solving the Einstein equations ^ in vacuum: 



G/ + 5/Ao = , 



(9) 



where Aq is given by ([5]) and ^ with T^^ = 0. For 
example, if f{R) = R — p'^/R, one finds Aq — ±\/3p'^/4:. 
It is well known that the solution to equation ^ is the 
static Schwarzschild-dS metric (or Schwarzschild-AdS if 
Ao < 0): 



Git) 1 



2GM A 



0^2 



= I 1 



2GAf Ao 



(10) 
(11) 



Here the integration constant M is time-independent 
since Gq^ = implies that dB/dt = Note that 

introducing the Lagrangian parameter G in pUj) and (fTT|) 
merely defines a convenient choice of units for M. The 
other integration constant G (t) is usually absorbed into a 
redifinition of t, but since we want to match the vacuum 
solution to a stellar interior, we will keep it explicit for 
the moment. 

The gravitational field is sourced by the trace of the 
energy-momentum tensor T = —p + 3p via ([5]). For all 
ordinary, non-relativistic stars, p <ti p (for example, at 
the center of the Sun p ~ 10~^p) and we can hence ne- 
glect pressure when solving for the metric. Moreover, we 
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are here considering an ideahzed object with a constant 
density, so that T — const. For this case the Einstein 
equations are again of the form ^ : 



(12) 



It is straightforward to show that the sohition to (fT2|) 
which is flat and non-singular for an observer at the origin 
is 



^ / 2Gm(r) Ap^^V'^' 



^ _ 2Gm{r) _ A,^2 
r 3 



where 



m(r) — d 



r 



(13) 
(14) 

(15) 



and Ap can again be found via ^ and ([7]) after the func- 
tion f{R) is specified. This completely fixes the stellar 
interior so that the only unknown quantities in the full 
solution are the parameters M and C (t) of the exterior 
vacuum. Matching at r = gives 



M = m(ro) -I- 



Ap - Aq ^ 
6G ' 



(16) 



8jrGp/Fp+Ap 



and C{t) = (l - 2GM/rQ - Aor|/3) 
i.e. a constant independent of time, which can be ab- 
sorbed by a simple scaling of t so that the exterior solu- 
tion is the standard Schwarzschild-dS metric. 

The possibility of finding a matching solution is of 
course obvious in the Palatini formalism. It follows di- 
rectly from the uniqueness of the exterior solution and 
the continuous nature of fixing i? as a function of T, 
both properties following directly from the trace equa- 
tion ([5]). However, what is not trivial is the nonstandard 
relation between the exterior mass parameter M and the 
interior density p, as shown in (|16p together with psp . 
Since the local density is what defines the local pres- 
sure and other thermodynamical properties of the star, 
it is obvious that F can not differ significantly from 1 
inside the Sun without changing the predictions of So- 
lar physics. Note that although the contribution to M 
due to the interior effective cosmological constant Ap will 
vanish in the limit -Fp — s- 1, the term with Aq will remain 
so that M will still differ from the GR value. This is 
simply due to the fact that F = 1 corresponds to GR 
without a cosmological constant, whereas in general the 
vacuum value of F = const. ^ 1 so that the theory cor- 
responds to GR with a cosmological constant Aq. Given 
that Aq correspond to the observed amount of dark en- 
ergy 17 A « 0.72 this subtle shift in mass is neverthe- 
less completely neglible. 

Interestingly, for all reasonable models for which 
SF{R) = F{R) — 1 decreases as a function of R, Fp tends 



very strongly towards unity. Indeed, for the particular 
case f{R) = R — n'^/R, we find from the trace equa- 
tion © that 



R=-[ 8nGp ± y/{8TTGp)^ + 12^4 



(17) 



in the Newtonian limit. If the exterior solution is taken 
to be asymptotically De Sitter space we must pick the 
solution with a positive sign in (|17p . Hence, given that 
fi'^ is scaled to cause the observed acceleration, then 
/i^ < SttGp and hence F - I = / R^ < 1 inside the 
Sun. In this case the changes to the stellar physics are 
completely neglible in a 1/i? model. Purely on academic 
interest, we note that the situation would be markedly 
different if one tried to match the f{R) interior solution 
with an external anti-De Sitter solution. Here, to match 
R smoothly to the limiting AdS-value, one would have 
to take the negative sign in ([T7|) . Then R would tend to 
zero and F become very large inside the star. This would 
imply that the maximum value for the exterior mass pa- 
rameter of a star with radius tq , is constrained to be on 
the order of 2GM /i^r|. 



III. THE TOLMAN-OPPENHEIMER-VOLKOFF 
EQUATION 

Let us now consider the general case of a relativis- 
tic star in static equilibrium with an equation of state 
P ~ p{p)- We start with the conservation law for energy- 
momentum which is of the regular form also in the Pala- 
tini formahsm of J{R) gravity 15]. That is, V ^T^"' = 0, 

where Vp is the covariant derivative with respect to the 
metric connection { p^ } . Assuming spherical symmetry 
and that we can describe matter inside the star by a per- 
fect fluid, T^"' = (p + p)uf'u'' - pgf"" , then gives 



A' 

V = -y(p + p) 



(18) 



where a prime denotes differentiation with respect to r. 
In GR one can use the regular Einstein equations to ex- 
press A' in terms of p, p and a mass function m(r) sim- 
ilar to (fT5| . but with Fp = 1, after which ([18]) gives the 
usual Tolman-Oppenheimer-Volkoff (TOV) equation for 
hydrostatic equilibrium. 

In order to find the generalized TOV equation for /(i?) 
gravity, we also need to solve A from the field equations 
(131111). From dH) it follows that 



i?M-(r) 



1 

F 



1 



VpV, - -.gp.D F , (19) 



where a tilde denotes that a quantity is defined in terms 
of the metric connection {^i^}. We can hence rewrite the 
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field equations ([3]) as 
inG. 



Giiu — 



3 1 



R- 



F 



- ( (V^F)(V.F)--g^.(Vi^)' 



F 



(20) 



wliere G^^ = i?^^ — \Rg^i,. Note tiiat / and F are still 
functions of the Ricci scalar R = g'^'^Rfj_^{r), and hence 
algebraic functions of T via the trace equation ([5]) . 

Using the 00 and 11 components of the field equa- 
tions (j20[) for a static, spherically symmetric metric, one 
obtains the following source equations for A{r) and B{r): 



A' 
B' 
where 



-1 
1+7 
1 

1 + 7 



1-e^ 



a 

—SirGrp + - 
r r 



e 



B 



-SnGrp + 



a 



(21) 
(22) 



a = r \ - 



— „2 




and 7 = rF' /2F. The TOV equation for f{R) gravity in 
the Palatini formalism now follows from (fTSl) and (1211): 



1 



(P + P) 



1 + 7 r(r - 2Gmtot) 

^ AirGr^p 
trmtot H 



(25) 



-(r - 2GTOtot) 



where we defined a new mass parameter 
2Gmtot(r-) = r(l-e-^). 



(26) 



Unlike in GR, it is in general not possible to obtain an 
explicit solution for B(r), or equivalently, mtotif) can no 
longer be found via a simple integration. However, we 
can use (|22p to obtain a differential equation for mtot (f ) : 



1 



F ' 
mtot 



r 



2G 



(27) 



where B is completely eliminated in favour of mtot (by 
eventually using ([26|) in the expression ([23|l for a) . Given 
the equation of state p{p), the TOV and mass equa- 
tion ((25l) and ((27|) will completely determine p(7'), /o(r), 
and mtot ('') j and hence also the metric via definition (I26p 
and the source equation (|2ip . The exterior mass parame- 
ter M is then given by Af = mtot('''0) — AoTq/GG, where 




2 3 4 

r/2GMQ 



xlO" 



FIG. 1: The effective cosmological constant A(r) (solid) and 
F{r) (dashed) for f{R) — R — fJ^'^/R and an assumed den- 
sity profile p(r) — po/(l + e^''"''©') (dotted), where po have 
been chosen so that M = Mq, the radius tq = 2.3 x 10* 
(3 km/2GM0), and ^ = 1.4 x lO"'^. Above, A and p have 
been normalized to the observed amount of dark energy 
flA ~ 0.72 [3| and the central density p(0), respectively. 



Aq is the background, vacuum value of the effective cos- 
mological constant. 

The true equation of state for matter inside a star has 
to be determined from micro-physical properties, and it 
is an external input to the gravity sector of the theory. 
We will not attempt to solve such a problem in this pa- 
per. Since pressure is neglible when computing the met- 
ric components of an ordinary, non-relativistic star, it is 
enough for the problem at hand to guess some smooth 
density profile p{r) and then directly use equations (PT|) 
and ((22)) to determine the metric [1^. For the particu- 
lar case of f{R) — R — p.'^/R the trace equation (O still 
implies that F « 1 inside the star, so that the stellar 
interior should look very much like in GR. One can see 
how this limit is approached from Fig. [1] where we plot 
F and the effective cosmological constant 



x^Ur-1 

2 V F 



(28) 



as functions of r for a p{r) that roughly represents the 
density profile of a star. Clearly, F = 1 inside the star 
and the contribution to AI coming from the effective cos- 
mological constant A will indeed be neglible. Obviously 
this example confirms the result of the previous section, 
that a 1/i? model is completely indistinguishable from 
GR at the Solar System scale. 



IV. SUMMARY AND DISCUSSION 

We have studied the interior spacetime of stars in the 
Palatini formalism of f{R) gravity. In particular, we de- 
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rived the generalized Tolman-Oppenheimer-VolkofF equa- 
tions for spherical hydrostatic equlibrium. In the Palatini 
formalism, the trace of the field equations provide an al- 
gebraic relation between R and T that gives a unique 
vacuum solution corresponding to the Schwarzschild-De 
Sitter metric 19]. Hence, in contrast to metric f{R) 
gravity, the interior spacetime of a star will not affect 
the form of the vacuum solution. This is the very reason 
for why Palatini f{R) gravity can pass the Solar Sys- 
tem constraints. However, matching the stellar interior 
with the exterior vacuum will give a non-standard rela- 
tion between the gravitational mass and the density of a 
star. For the particular case /(i?) = R — l^'^/R, the ef- 
fective cosmological constant A vanishes inside the star, 
i.e. F ^ 1, so that the addition of a 1/R term to the 
Einstein-Hilbert action is consistent with Solar System 
constraints. This is due to the fact that the curvature R 
becomes large compared to /i^ inside the star, suppress- 
ing the effect of the fi'^/R term. However, if one considers 
for example a term eR^, its effect would grow inside the 
star. One might thus see interesting modifications to the 
physics of compact stars, and thereby perhaps be able to 



put constraints on parameters such as e. 

In summary, although one can arrange for different 
forms of f{R) to give rise to the same cosmological con- 
stant in the exterior vacuum solution, the corresponding 
stellar interior solutions can be vastly different depend- 
ing on the form of f{R). This raises the possibly to 
differentiate between different f{R) models and General 
Relativity, depending on their effect for the stellar evo- 
lution and stellar paramers such as luminosity. However, 
neither Solar System tests nor astrophysical considera- 
tions are able to rule out the addition of a 1/R term to 
the Einstein-Hilbert action for gravity. 
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